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We study the few-electron eigenspectrum of a nanotube quantum dot with spin-orbit coupling. 
The two-electron phase diagram as a function of the length of the dot and the applied parallel 
magnetic field shows clear signatures of both spin-orbit coupling and electron-electron interaction. 
Below a certain critical length, ground state transitions are correctly predicted by a single-particle 
picture and are mainly independent of the length of the dot despite the presence of strong correla- 
tions. However, for longer quantum dots the critical magnetic field strongly decreases with increasing 
length, which is a pure interaction effect. In fact, the new ground state is spin- and valley-polarized, 
which implies a strong occupation of higher longitudinal modes. 



I. INTRODUCTION 

Carbon nanotubes have allowed to realize clean quasi- 
one-dimensional electron systems. Experiments reveal- 
ing fundamental interaction effects include the detec- 
tion of Wigner crystallization^ or Luttinger liquid like 
behavior^''^. An interesting feature of nanotubes is that 
the orbital part of low lying excitations has an additional 
spin-like degree of freedom the valley index. This new de- 
gree of freedom can cause orbital Kondo effect** unusual 
spin configurations^, or a new type of shell structure^ in 
nanotube quantum dots. 

It was generally assumed that spin and valley de- 
grees of freedom lead to a fourfold degeneracy of elec- 
tronic states, however, recently spin-orbit coupling was 
observed to split this degeneracy in two pairs of either 
parallel and antiparallel spin and valley orientation^. In- 
terestingly the experimental data could be well explained 
in a single-particle picture, and correlation effects seemed 
to be of minor importance. In this work we analyze how 
interaction effects show up in the two-particle spectrum 
of a single nanotube quantum dot with spin-orbit cou- 
pling. We argue that the eigenspectrum can be divided 
in multiplets of states that have the same orbital symme- 
try. Energy gaps within the same multiplet are only de- 
termined by spin-orbit coupling and the orbital Zeeman 
effect (and additional small correction due to local inter- 
actions), and are therefore captured in a single-particle 
picture. However, the extend of correlations can be ap- 
preciated by comparing different multiplets. In partic- 
ular we show that above a certain critical length a tiny 
magnetic field is enough to cause a ground state tran- 
sition to a spin and valley-polarized two-particle state, 
that necessarily involves the occupation of higher modes. 

In the next section we introduce our model. The 
quantum dot is described by a potential well along the 
nanotube and a continuum description is applied for the 
single-particle spectrum of electrons localized in this well 
and subject to a parallel magnetic field^. The single- 
particle spectrum also includes the effect of spin-orbit 
coupling^'^. We then show how the electron-electron in- 
teraction can be correctly incorporated in the continuum 
model^'*°. Thereafter we present our results, including 



a detailed discussion of the phase diagram of the two- 
electron ground state as a function of magnetic field and 
length of the quantum dot. 



II. MODEL 

A nanotube is a one-atom-thick layer of graphite called 
graphene wrapped into a seamless cylinder. Depending 
on the orientation of the underlying honeycomb lattice 
of carbon atoms with respect to the symmetry axis of 
the nanotube, it is either metallic or semiconducting*^. 
We will study a semiconducting nanotube with an ad- 
ditional confinement potential along the tube, which is 
controlled by external gates and gives rise to a discrete 
set of localized electronic states. 



A. Single-particle spectrum 

In a continuum description, the single particle orbitals 
have two components belonging to the two sublattices, 
called A and B in the following. Furthermore, the single- 
particle states have an additional spin-like degree of free- 
dom T S ±1, the valley index, since there are two in- 
equivalent band minima at the K and K' = — K points 
of the graphene's Brillouin zone. 

Using cylindrical coordinates C, <j) the single-particle 
Hamiltonian is given by 

Ho = -ihvFira^^d^ + aydc) + ViO , (1) 

where vp is the Fermi velocity and ax , (Jy are Pauli ma- 
trices acting on the sublattice space. We study a square 
well potential, i.e. V{Q is zero for \(\ < L/2 and Vq 
otherwise*. We assume the potential to be smooth on the 
atomic length scale (interatomic distance aq — a/y/3 = 
0.142 nm, where a is the lattice spacing) and therefore 
neglect confinement induced intervalley scattering. 
The single particle solutions are given by: 



^r.kir) = (2^i?)-i/2e'^K'-e^««^0fe(C) 



(2) 



where fc, k denote the wavevectors along and around the 
tube and the two component longitudinal wavefunction 
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is normalized such that J d( {\(j)Ak\'^{C) + |0BfcP(C)) = 
1. </'fc(C) is given by a standing wave with wavevec- 
tor k inside the weU and evanescent modes outside the 
weh*. The corresponding eigenenergy is given by Ek = 
hvF\/ + hP' . We note that we measure energy with 
respect to the center of the gap, so that the dominant 
part of the single-particle energy is constant and given 
by hvpK ~ 220meV/ R[nm\. Electron-electron interac- 
tion however, affect the longitudinal part with a much 
smaller level spacing that depending on the length of the 
dot is A = 2 - 10 meV. 

Both the axial magnetic field B and the spin-orbit cou- 
pling modify the transverse wavcvector k 



Ago = 0.37 meV 
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K t/R{1/'3 + T$/$o + rfT$5o/$o) 



(3) 



Here a denotes the spin component along the tube, $ = 
ttR^B the magnetic flux through the tube, $o = h/e, 
and <&so ~ 7.2 10~^ determines the curvature induced 
spin-orbit interaction^'^. The second term in Eq. (3) re- 
sults from the coupling between the external magnetic 
field and the orbital magnetic moment that is caused by 
the transverse motion around the tube^^. Electrons in 
different valleys have an opposite sign of this orbital mo- 
mentum, which leads to a valley splitting that is linear 
in the applied magnetic field. We call this the orbital 
Zceman term in analogy with the smaller spin Zeeman 
term Hz = —q^bctB /2 that leads to spin dependent 
energy shift in the magnetic field. The orbital mag- 
netic moment is /iorb ~ — T0.5i?[nm] meV/T and the 
spin magnetic moment ^spin ~ (t0.06 meV/T . The 
third term in Eq. (3) describes spin-orbit coupling. It 
increases (decreases) the energy of single-particle states 
with aligned (anti-aligned) spin and valley degree of free- 
dom. The energy splitting is approximately given by 
Aso ~ 1.9/(i[nm] meV, where d = 2R denotes the di- 
ameter of the tube.*" 

Spin-orbit coupling and orbital Zeeman effect couple 
to the transverse part of the wavefunction while their ef- 
fect on the longitudinal part (j)k{0 can be neglected for 
the large dot sizes we are interested in. The longitu- 
dinal wavevector k is determined by the transcendental 
equation^ 



tan(fcL) 



kk 



Ek{Ek - VG)l{hvpY - 



1/2 



(4) 



where k — — \^E\^ — Vg) I^'^fX ^ determines the 
decay of the wavefunction outside the well. Due to the 
symmetries of the Hamiltonian in Eq. (1), (t>k{C) is real 
and has a well-defined parity p — ±1, 0a(C) = p4'b{—C)^ 
where A, B label the two sublattices. The parity of the 
j-th mode (where the ground state corresponds to i = 0) 
is given hy p — (— !)'■ 

Figure 1 shows the magnetic field dependence of the 
two lowest longitudinal modes. Each mode gives rise to 
four single-particle states due to the two spin and two 
valley degrees of freedom. At zero magnetic field these 
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FIG. 1: (Color online) Single-particle states for a, R = 2.5nm 
nanotube with a. L — 80nm square well of depth Vg = 50meV . 
Color coding: Red (solid): t = —, a —1, Green (long dashed): 
T = +1 a =1, Blue (short dashed): t = —, a Purple 
(dotted): t — +, a —'\. 



four states are split in two Kramer doublets (states ob- 
tained by flipping simultaneously spin and valley degree 
of freedom are degenerate due to time-reversal symme- 
try). Finite magnetic fields lead to energy shifts linear in 
magnetic field caused by orbital and spin Zeeman split- 
ting. Denoting the single-particle states by P, r, a, a 
single-particle picture predicts the two-particle ground 
state to be |-I-1,-,T; for B < B^^t « 0.15T 

and 1+1,— ,t; +1,— ,J,) for B > Bcrit- We note that 
only for ridiculously large magnetic fields of about 80 
Tesla the single-particle picture predicts a further ground 
state crossing due to an occupation of the first excited 
shdl |-H,-,T; -1,-,T)- 



Interaction 



Due to the large energy gap between different trans- 
verse modes, it is justified to treat completely filled as 
well as completely empty transverse subbands as inert, 
giving rise to a static screening constant e. Assuming 
gate electrodes to be sufficiently far away from the quan- 
tum dot, we use a long-ranged interaction between the 
conduction electrons ?7/(ri,r2) = e^/(e|ri — r2|). 

The interaction does not depend on the electron spin 
and is therefore diagonal in the spin degree of freedom. 
However, the local part of the interaction is not diago- 
nal in the valley degree of freedom'^' It is therefore 
instructive to split the interaction Hint in a long-ranged 
part Vc and a local, onsite interaction Vh- 

After integrating out the transverse motion the inter- 
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action is given by: 

Hint ~ Vc +Vh 

= ^ 1 dCi dC2V{Ci - C2)P(C1)P(C2) 



a) AgQ = 0.37 meV, a = 



b) AgQ = 0.37 meV, a = 1 



Here p(C) = Ep^r V'ji<Tr(C)'0p<Tr(C) labels the charge 
density at (, where p G {A, B} denotes the sublat- 
tice index. The field operators are now expressed by 
^JpariC) — J2k (t^pkiOdark whcrc ttark annihilates an elec- 
tron in the longitudinal mode k (characterized by the 
wavefunction (j)k of Eq. (2)) and with spin a and val- 
ley T. The now one-dimensional interaction V{C) = 
2e2X [4i? V(C2 -|-4i?2 )] / [g7r(,^2 _^4i?2 ) i/2] can be expressed 
by the incomplete elliptical integral of first kind K(x)^'^. 
In the following the strength of the long-ranged interac- 
tion is characterized by the effective fine structure con- 
stant a — /ehvF ~ 2.2/e. The local part depends on 
U = U Au/ {2t:R), where Au = a^-\/3/2 denotes the size of 
the graphene's unit cell in real space and U is the onsite 
interaction. We use U = 15eV.^ 

Due to the rapidly oscillating Bloch factors e*"^^"" of 
the eigenfunctions (2), the long-ranged interaction Vc is 
diagonal in the valley and spin degrees of freedom, and 
in agreement with the continuum description, the inter- 
atomic distance between the two sublattices is neglected. 
The lattice effects not captured in the continuum model 
and the long-ranged interaction Vc are taken into ac- 
count by the local part of the interaction Vh- We note 
that Vc depends equally on spin- and valley symmetry, 
but this is not the case for local interaction. For ex- 
ample spin-aligned electrons do not interact via V//, but 
valley-aligned do. While local interactions can be very 
important for short quantum dots^"* their effect is rather 
small for the long quantum dots as shown in the fol- 
lowing. However, also the spin-orbit coupling is a small 
quantity and as we will discuss below local interaction en- 
ergies can add up to the spin-orbit interaction. We note 
that Vh still conserves valley polarization and that 

it allows for intervalley exchange interaction (V^ ' ) . 

The many body eigenfunctions can be characterized 
by a triple of quantum numbers (P,Sz,Tz), where P ~ 
YlnPn ^ {=1=1} denotes the total parity, Sz — l/2X]n^n 
the z- component of the total spin and Tz — 1/2 X^n'''" 
the total valley polarization; where n = 1, .., runs over 
all electrons. We note that without local and spin-orbit 
interaction the two-particle states can also be chosen as 
eigenstates of total spin S*^ and total valley degree of 
freedom T^. 

In the following we calculate the few-electron eigen- 
spectrum of the Hamiltonian H = Hq + Hz + Hint ■ We 
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FIG. 2: (Color online) Energies of lowest two-particle states. 
Left no interaction, right interaction a = /ehvp ~ 1. Pa- 
rameters as in Fig 1. Blue (solid): Parity -1-1, Red (dashed): 
Parity —1. Three states are marked: Crosses: (P = 1,T^ = 
0, Sz =0), Filled squares: (1, —1, 0), Open circles: (1, —1, 1). 



restrict the single-particle basis to the bound longitudi- 
nal modes of the lowest transverse mode and diagonalize 
the few-electron Hamiltonian for each set of conserved 
quantum numbers {P,Tz,Sz), so that electron-electron 
correlations within this basis set are fully taken into ac- 
count. 



III. RESULTS 

Figure 2 shows the two-particle spectrum as a func- 
tion of magnetic field both for non-interacting (left part 
a — 0) and interacting electrons (right part a = 1). For 
non-interacting electrons the ground state corresponds to 
a double-occupation of the lowest longitudinal mode and 
has parity P = -1-1 (blue states). Since the spacing A to 
the next longitudinal mode is much larger than the spin- 
orbit splitting, states with negative parity P = — 1 (red 
states) are energetically well separated from the ground 
state for all relevant magnetic fields. In the following we 
label eigenstates by their quantum numbers {P,Tz, Sz)- 
Without magnetic field the nondegenerate ground state 
corresponds to the subspace with (1, 0, 0) (blue line with 
crosses) which is favored by spin-orbit coupling. At a crit- 
ical magnetic field the ground state crosses to (1,-1,0) 
(blue line with filled squares) due to the orbital Zeeman 
term. 

Electron-electron interaction strongly reduces the gap 
between the P = 1 and P = — 1 states as shown on the 
right side of FIG. 2. At the same time spin-orbit induced 
energy gaps are unaffected by interactions and also the 
magnetic field dependence of the energies is the same as 
in the noninteracting case. For the parameters chosen 
in FIG. 2 the ground state transition at finite magnetic 
field occurs to the (—1,-1,1) state, which is spin and 
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b) U = 15eV, Aso = 0.37meV 



T 



T 



) l ( I ) l ( ) l ( ± ) l ( )K I ) l ( )T( 



0.4 



0.6 



0.8 



FIG. 3: (Color online) Dependence of two-electron spectra on 
strength of long-ranged interaction a in absence of a mag- 
netic field for R — 2.5nm, L — lOOnm. a) Total two-particle 
energies in absence of local interactions and spin-orbit cou- 
pling. Colors label symmetry of the longitudinal part of the 
wavefunction. Blue (solid) and green (short dashed) indi- 
cate symmetric and red (long dashed) an antisymmetric lon- 
gitudinal wavefunction. b) Two-particle excitation energy 
AE = Ei{N ^ 2) - Eo{N = 2). Energy splitting within 
each multiplet is determined by spin-orbit coupling and local 
interactions. Meaning of marked states as in FIG. 2. 



valley-polarized (red line with open circles in FIG. 2). 

We now study the reduction of the energy spacing 
between the P = -1-1 and P = —1 multiplet of states 
with increasing interactions. Since the magnetic field de- 
pendence of the energies is hardly changed by interac- 
tions it is instructive to study the spectrum in absence 
of magnetic fields. We first neglect spin-orbit coupling, 
^SO = and the onsite interaction, U = 0. Figure 3 a) 
shows the eigenspectrum oi Hq + Vc a.t B — as func- 
tion of the interaction strength a. The blue (solid) line 
is the 6-fold degenerate P = +1 ground state. Inter- 
action split the P = —1 states in two sets. The red 
(long dashed) line consists of 10 states that approach 
the P = +1 ground state with increasing interactions. 
The green (short dashed) line indicates the remaining 
six states with P = —1. Vc exclusively acts on the longi- 
tudinal part of the wavefunction since a single transverse 
subband is considered. The two-particle eigenstates can 
be factorized in longitudinal, spin and valley parts. The 
longitudinal part of the P = +1 ground state is sym- 
metric with respect to interchange of two electrons. It 
is multiplied with either valley-triplet and spin-singlet 
or valley-singlet and spin-triplet, in order to guarantee 



the antisymmetry of the total two-particle wavefunction. 
This explains the six-fold degeneracy. The energetically 
favoured set of P = — 1 states is ten- fold degenerate and 
has an antisymmetric orbital part and spin and valley 
part are either both singlet or both triplet. Figure 3 a) 
shows that the two-particle ground state has always a 
symmetric orbital part for all interaction strengths, how- 
ever the energy gap to the eigenstates with antisymmet- 
ric longitudinal part vanishes with increasing interaction. 
We note that for scalar eigenfunctions of a Schrodinger 
equation a symmetric orbital part is guaranteed by the 
Lieb-Mattis theorem^^. Since we are describing a semi- 
conducting nanotube with a large gap between transverse 
modes, we are in fact very close to that limit. With in- 
creasing a, electrons become more and more correlated 
and finally form a quasi classical Wigner crystal, where 
electrons are localized in different region of space and 
symmetry becomes irrelevant. 

Including again spin-orbit coupling and local inter- 
actions states with an originally symmetric (antisym- 
metric) longitudinal part split in multiplets of six (ten) 
states. This splitting is small with respect to the total 
energy, so that Figure 3 b) shows the excitation energy 
AE = Ei{N = 2) - Eo{N = 2) rather than the total 
energies. The two-electron ground state, which belongs 
to (-f 1,0,0), therefore defines AE = axis. Excitations 
toP = -|-l(P = — 1) states are depicted in blue (red). 
Since an increase of the interaction strength a leads to an 
increasing distance between the two electrons the proba- 
bility of finding both on the same site strongly decreases 
with increasing a and the local interaction effects vanish. 
The energy splitting within the different multiplets there- 
fore quickly approaches the constant spin-orbit gap with 
increasing a. Effects of local interactions are however vis- 
ible if the quantum dot becomes shorter, or if a is small 
(we assume that local interactions are not screened) . On- 
site interactions favor spin triplet states over spin singlet 
states and increase the energy gap between the (1,0,0) 
ground state, that is in a superposition of spin singlet 
and spin triplet, and the (1, —1, 0) state (filled squares in 
FIG 3b)), which is a spin singlet state. 

Above a critical magnetic field, a valley-polarized 
ground state (T^ = —1) is favored, the parity of which 
depends on the ratio of single-particle and interaction 
energy. This ratio increases for decreasing length of the 
quantum dot or increasing radius of the nanotube or de- 
creasing dielectric constant e. The length of the quantum 
dot is tunable experimentally by changing gate voltages. 

We now discuss the phase diagram of the two-electron 
quantum dot as a function of magnetic field P, and length 
L of the quantum dot, for different spin-orbit couplings 
Aso and interaction strengths a as shown in FIG. 4. We 
note that the appearance of the spin and valley-polarized 
state (green areas) is favored by the interplay between 
spin-orbit coupling and long-ranged Coulomb interac- 
tion. Figure 4 a) and b) show the phase diagram without 
spin-orbit coupling. Then the ground state in zero field 
is given by the three spin-polarized states (among which 
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a) AsQ=0, a = 0.3 b) Asq=0, a = 1 . 
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FIG. 4: (Color online) Two-particle ground state as a function 
of magnetic field and length of the nanotube for R = 2.5nm. 
The ground state belongs to one of the following sets of 
quantum numbers: Blue: (1,0,0) Red: (1,-1,0), Green: 
(— 1, —1, 1). Upper (lower) row shows phase diagrams in ab- 
sence (presence) of spin-orbit coupling and left (right) column 
shows different interaction strengths of long-ranged Coulomb 
interaction. 
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FIG. 5: (Color online) Dependence of critical length L on 
interaction strength a and radius R of the tube. If the quan- 
tum dot is longer than the critical length the spin and valley- 
polarized state (—1, —1, 1) becomes ground state at a small 
magnetic field. Spin-orbit coupling Aso = 1.9 meV/d[nm]. 



a) Ago = 0.37 meV, L = 60nm b) AgQ = 0.37 meV, L = 1 0Onm 
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FIG. 6: (Color online) Energy needed to add a second particle 
to the quantum dot. Only transitions that are allowed by spin 
and valley selection rules are shown. R = 2.5nm, a — 1. 



the (1, 0, 0) one is indicated by the blue area on Fig 4 ). 
They are separated from the three spin singlet states of 
the P — 1 multiplet by local interaction effects only. At a 
critical magnetic field the valley-polarized state (1,-1,0) 
(red area of Fig 4) is favored due to the orbital Zee- 
man term. In agreement with our discussion above, the 
local interaction is more relevant for shorter quantum 
dots and for small a. For sufficiently large quantum 
dots the (—1, —1, 1) state (green area of Fig 4) becomes 
the ground state since long-ranged Coulomb interaction 
strongly suppresses the level spacing between the P = 1 
and P = —1 states and the remaining gap can be com- 
pensated by the gain in the (spin) Zeeman term. 

Figure 4 c) and d) show that the phase diagram dras- 
tically changes if spin-orbit coupling is included. The 
zero-field ground state (still belonging to (1,0,0)) is now 
non-degenerate and is in a superposition of spin singlets 
and triplets. Additionally, the regions where the ground 
state belongs to either (1,0,0) or (—1,-1,1) are both 
considerably enlarged at the expense of the (1,-1,0) 
state. Below a critical length, the magnetic field where 
the ground state crossing (from (1, 0, 0) to (1, —1, 0)) oc- 
curs is mostly length independent and coincides with the 
value predicted by a single-particle picture^. In con- 
trast for somewhat larger quantum dots the transition 
occurs to the (-1,-1,1) state and the corresponding crit- 
ical magnetic fields decrease continuously with increasing 
length of the quantum dot. 

In the phase diagrams of the two-particle ground state 
with spin-orbit coupling there is generally a critical 
length above which a magnetic field causes a ground state 
transition to the spin and valley-polarized state. The de- 
pendence of this critical length on a and the radius of 
the nanotube is shown in FIG. 5. In agreement with our 
discussion this critical length decreases with increasing 
interaction strength or decreasing radius. 

A powerful tool to measure the few-electron spectrum 
of the quantum dot is transport spectroscopy for differ- 
ent lengths of the quantum dot. Such an experiment al- 
lows to measure the energy needed to cause a transition 
from the one-electron ground state of energy Eo{N = 1) 
to a two-particle excited state Ei{N = 2) where i de- 
notes the excitation. Allowed transitions from the one- 
particle ground state to a two-particle excited state can- 
not change Tz or Sz by more than ±1/2. These tran- 
sitions are depicted in FIG. 6. For the shorter quan- 
tum dot, the energy needed for the iV = 1 to = 2 
ground state transition exactly follows the first excited 
one-electron energy (except for a constant charging en- 
ergy). This is not the case for the longer quantum dot, 
where the two-electron ground state transition occurs at 
smaller fields. We note that excitations between two- 
particle states of different parity are always modified by 
interactions and are not a mere combination of level spac- 
ing and spin-orbit gaps. 
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IV. CONCLUSIONS 

We have presented a detailed study of the two-electron 
cigcnspcctrum of a nanotubc quantum dot with spin- 
orbit coupling. Generally we find that the eigentates are 
strongly correlated and by varying the length of the quan- 
tum dot we identify clear signatures of short and long- 
ranged interaction. In particular we studied the two- 
electron phase diagram as a function of the length of the 
quantum dot and the applied magnetic field. While the 
ground state at zero magnetic field always corresponds 
to the same set of quantum numbers (given by parity 
P = 1, spin Sz = and valley polarization = 0) 
for all lengths, a finite magnetic field causes a transi- 
tion to a valley-polarized state which is cither a spin sin- 
glet or a spin triplet, depending on the length of the 
quantum dot. The former case is the one predicted by a 
single-particle picture since valley-polarized electrons in 
the lowest mode must be in a spin singlet state. This 
crossing is unaltered by Coulomb interaction, since the 
two crossing states have the same orbital part and their 
splitting is only given by the length independent orbital 
Zeeman shift and the spin-orbit gap (plus small correc- 
tion due to local interactions). However, once the length 
of the quantum dot exceeds a certain critical value inter- 
action effects cause a ground state transition to a spin 
and valley-polarized state. Increasing the length even 
further, the magnetic field of the ground state transition 



becomes arbitrarily small. 

An interesting continuation of our work is to analyze 
interaction effects on recently suggested optical manipu- 
lation schemes for the spin in carbon nanotubes^^. An- 
other open question is whether the presented interaction 
effects arc also manifest in double quantum dots in nan- 
otubes. In particular we suggest to study whether a val- 
ley and spin polarized ground state modifies selection 
rules leading to spin and/or valley blockade in transport 
spectroscopy^^. 

Note added While preparing this manuscript we be- 
came aware of the work of A. Secchi and M. Rontani [18], 
who obtained similar results for a nanotube quantum dot 
with harmonic confinement instead of the potential well 
used here. We note that the agreement of general conclu- 
sions in both works shows the robustness of the discussed 
interaction effects. 
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